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Q\ ' Abstract 

O ' 



We describe the statistical behavior of anisotropic crystalline membranes. In 
particular we give the phase diagram and critical exponents for phantom mem- 
branes and discuss the generalization to self-avoiding membranes. 



Jj^ 1 1 Introduction 



The statistical mechanics of phantom tethered membranes with bending rigidity is quite 
^ ■ well understood if the elastic and bending moduli are isotropic. There is a continuous crum- 

pling transition from the expected high-temperature crumpled phase to a low-temperature 
flat phase that has no analogue in conventional spin systems |l|, |2| . This flat phase is char- 
acterized by long-range orientational order of the membrane in the embedding space and 
large longitudinal (in-plane) fluctuations. Two important generalizations of this class of 
membranes arise by incorporating self- avoidance and/or intrinsic anisotropy. In ||] Radz- 
ihovsky and Toner (RT) showed that intrinsic anisotropy generally leads to a much richer 
phase diagram for phantom (non-self-avoiding) tethered membranes. Although anisotropy 
is irrelevant in the flat and crumpled phases it leads to an entirely new tubular phase sep- 
arated by distinct continuous phase transitions from both the conventional phases. This 
phase diagram is shown in Fig. p]. 

^From the applied point of view both isotropic and anisotropic membranes should be 
important. Isotropic membranes may be made by suitable random polymerization of fluid 

1 Research supported by the Department of Energy U.S. A under Contract No. DE-FG02-85ER40237, by 
Syracuse University by Alexander von Humboldt Stiftung, and by the Deutsche Forschungsgemeinschaft. 
2 Talk given at Lattice 98, Boulder Colorado, July 1998. 



1 




Figure 1: The mean field theory phase diagram for phantom anisotropic membranes from 



vesicles. On the other hand polymerization in the presence of an applied electric field 
should in principle yield an anisotropic membrane Q|. Polymerized membranes with in- 
plane tilt order have a natural anisotropy. The tubular phase is characterized by long-range 
orientational order in one (extended) direction only — in the remaining transverse direction 
the membrane is crumpled. 

The tubular phase thus resembles a rough sausage with a crumpled cross-section. This 
existence of this totally tubular phase was confirmed by Monte Carlo simulations in ||. 
In this paper the size (Flory) exponent up, describing the growth of the tubule-diameter 
Rq with internal system size (L), and the roughness exponent £, characterizing the height 
fluctuations along the extended tubule axis, were also reported. The results up = 0.305(14) 
and £ = 0.895(60) were in rough qualitative agreement with the theoretical predictions 
v = \ and £ = 1 of 0. More extensive subsequent simulations || have improved these 
first results - the best current results are up = 0.269(7) and Q = 0.859(40). Configurations 
characteristic of the various phases are shown in Fig. |2[ 



2 Model 



A tethered membrane is described by a 2-dimensional regular triangulated net, with the 
topology of a disk. Each node of the net has six neighbors, except at the boundary. The 
Hamiltonian of the system is 
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Figure 2: The three phases of anisotropic membranes: (a) tubular (b) crumpled and (c) flat. 
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where r CT is the position in 3-dimensional space of the node labelled a = (o~ x , a y ). The first 
sum runs over all nearest-neighbor pairs (bonds) of the membrane, and is the tethering 
potential. The second and third term are the bending energies in the x and y intrinsic 
directions. The bending energy is a ferromagnetic interaction between the unit normals to 
the faces of the membrane. The strength of this interaction is anisotropic: if two adjacent 
faces share a bond parallel to the x direction the coupling is n%; otherwise it is k-i- 
The canonical partition function for a membrane of fixed number of nodes N is 

Z = |[dr] 5(r cm ) e-^K (2) 

The Hamiltonian of Eq.(^) was simulated using Monte Carlo methods for triangular 
lattices ranging from 25 2 to 175 2 nodes. Further numerical details are given in [|| : H|. 



3 Global Phase Diagram 

To determine the global phase diagram of anisotropic membranes we explored the line 
(ki,K2) = (3k, k). We have also simulated the lines (k, 0) and (2, k) on smaller lattices. 
Distinct signatures of phase transitions were found in the specific heats Cy associated with 
the variance of the two bending energy terms in the action E x and E y ; 

CvM = p §~ K {Ei). (3) 

In Fig. H we show the total specific heat along the line k% = 3k2 for lattice sizes up to 175 2 . 
There are two distinct peaks. We have confirmed that the two peaks occur at different 
values of the bending rigidity, signaling the existence of two distinct transitions. 

Performing a similar analysis along the other two lines in the (^1,^2) plane yields a 
phase diagram consistent with that of Fig. [j]. This confirms the three phase scenario. 
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Figure 3: The value of the total specific heat Cy for various lattice sizes. The interpolating 
lines are obtained using multi-histogramming methods. 



4 Tubular Exponents 

A more detailed understanding of this tubular phase is obtained by looking at the fluctu- 
ations of the zero- mode of the tubule height h rms , analogous to the height fluctuations of 
a flat membrane, and the scaling of the width of the tubule R± |], 0. These are expected 
to scale as h rms ~ and R!]_ ~ L vp . 

We find £ = 0.859(40) and vp = 0.269(7). The result for up is in excellent agreement 
with the analytic continuum prediction of |j] {vp = 1/4) but the roughness exponent is 
well below the analytic prediction Q = 1. Further studies are under way Q. 

The most challenging problem in the physics of tubules is the incorporation of self- 
avoidance. Self-avoidance becomes relevant below the upper-critical embedding dimension 
d = 11, and here one can develop a systematic e expansion |K|. This expansion can be 
improved and this is treated in the contribution to this proceedings by Alex Travesset 
p| , |i~2| . We are also performing Monte Carlo simulations of this system. Preliminary 
results indicate a flat-to-tubular transition in the physical case of three dimensions. The 
specific heat plot is shown in Fig. |4|. 
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Specific heats for anisotropic self-avoiding membrane 
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Figure 4: The specific heat Cy for self-avoiding membranes. 
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